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Abstract 

Gross has found an exact expression for the density of eigenvalues in the simplest 
version of the Kazakov-Migdal model of induced QCD. In this paper we compute 
the spectrum of small fluctuations around Gross's semi-circular solution. By solving 
Migdal's wave equation we find a string-like spectrum which, in four dimensions, cor- 
responds to the infinite tower of mesons in strong coupling lattice QCD with adjoint 
matter. In one dimension our formula reproduces correctly the well known spectrum 
of the hermitean matrix model with a harmonic oscillator potential. We comment on 
the relevance of our results to the possibility of the model describing extended objects 
in more than one dimension. 



Introduction 



The original Kazakov-Migdal (KM) model of induced QCD [|T| consists of a scalar in the 
adjoint representation of SU{N), covariantly coupled to gauge fields on a hypercubic lattice 
of spacing a. There is no kinetic term for the gauge fields and the action of the model is 



s=ny: 



(1) 



The hope was that this model could perhaps induce QCD in the sense that a distances 
much larger than a a kinetic term for the gauge fields would be generated through scalar 
interactions at distances on the order of the cutoff a. By carefully tuning the parameters of 
the potential it might then be possible to reach an 'asymptotic freedom domain' in which 
the continuum limit could be taken. The mass of the scalar, which is kept heavy in this 
limit, would act as an effective cutoff for the resulting continuum (QCD?) theory. 

The beauty of this idea of course lies in the fact that the model given in Eq. (1) is 
analytically tractable in the limit of an infinite number of colors @. The model is 
especially simple in the case of a purely quadratic potential and initially there was clearly 
no a priori reason for why the objective of inducing QCD, if possible at all, could not be 
achieved by just tuning the mass of the scalar. This simplest of all possibilities is by now 
ruled out. The reason is as follows: D. Gross @ has found an exact solution to the saddle 
point equations describing the translationally invariant eigenvalue distribution of the $-field 
in the N = oo limit. This solution has no continuum limit in four dimensions. It rather 
describes the large A^, infinite coupling (remember, there is no plaquette term for the gauge 
fields in (1)) limit of lattice QCD with adjoint matter. This in by itself is not sufficient to 
prove 'non- induct ion' for there could be other saddle points. However computer simulations 
of the SU{2) theory |^ also show that there is no critical point in the case of a purely gaussian 
potential. Assuming that the large limit is smooth the theorem is proved. Incidentally, 
the simulations also showed that by adding a quartic term a critical point could be reached. 
The nature of the resulting continuum theory is not clear at the moment. Migdal |P has 
made considerable progress in the investigation of the theory with a general potential but 
lately he has concentrated his efforts on the study of an extended version of Eq. (1), the so 
called 'mixed model' of induced QCD 0. 

In the mixed model ^ Uf <^ N heavy fermions in the fundamental representation are 
added to the action in order to break the additional local symmetries of the model, in 
particular the local Zj^ symmetry U^{x) — > Z^{x)U^{x) first discussed by Kogan et al. [0. 
As pointed out by these authors, due to this symmetry the Wilson loop cannot acquire an 
expectation value. The adjoint loop of course does, but it is always screened and cannot 
serve as an order parameter for confinement. It is Migdal's hope that there is a phase of the 
mixed model in which the center symmetry is broken and Wilson loops show an area law. 

In the present paper we will have nothing more to say about the question of whether 
and how QCD can be induced. Rather we will extend the very pretty work of Ref. on 
the original Gaussian KM model by calculating the spectrum of small fluctuations around 
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the Gross saddle point. We believe that there are many reasons for why this is interesting. 
In four dimensions our result gives the masses of the mesons of the infinite coupling lattice 
gauge theory to leading order in The large A^, strong coupling spectrum of the theory 

with fundamental fermionic matter has been known for over a decade now and the adjoint 
scalar spectrum nicely complements this old work. Furthermore, the model in Eq. (1) can 
be viewed as a gauged matrix model. Hermitean matrix models in Z) < 1 dimensions are 
intimately related to discretized versions of the Polyakov string [Q] and as such describe 2D 
gravity coupled to matter fields. In particular, a case of great current interest is the one 
dimensional case, describing c = 1 matter coupled to 2D gravity. In more than one dimension 
ordinary matrix models become very complicated and one might hope that the gauged 
version in Eq. (1) might be useful as theory of extended objects in more than one dimension. 
It has also recently been sugg ested ^ that the KM model can be viewed in D dimensions 
as the high temperature limit of the {D + 1) dimensional Wilson action. If this intersting 
conjecture is indeed correct then our result in three dimensions should describe the spectrum 
of fluctuations of the 4D-Polyakov line around one of it's minima at infinite temperature. 
Finally, our calculation can serve as a testing ground for the methods developed by Migdal. 
In some sense the present problem represents the 'hydrogen atom' of induced QCD; we can 
actually work out the eigenvalues and the eigenf unctions of Migdal's wave equation in closed 
form. In the process we also found a term in the wave equation which had been omitted in 
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The Gross Solution 



In this section we will quickly review Gross's solution of induced QCD with a quadratic 
potential. In the process we will also establish notation. We will follow Migdal's approach 
rT[| to the problem and express the action in Eq. (1) in terms of the density of eigenvalues 

p(^) = ^E^(^-^^)- (2) 

i 

This is easily done by first integrating out the gauge fields using the Itzykson-Zuber integral 
which produces 

S = Nj2U{Xi{x)) - 5] ln|Ai(x) - Aj(x)| -^ln[/(A(x),A(x + H] (3) 

x,i x,iy^j x,fj. 

where the logarithm comes from two factors of the Vandermond determinant and I denotes 
the Itzykson-Zuber determinant. Now using Eq. (2) the action can be written in terms of a 
x-dependent density as 

S[p] = N^Yl fd\p^i\)Ui\)-N''J2 jdXpM f dX'p,{X')\n\\ - \'\-J2ln[Iip,,p^+^a)]. 

(4) 

Note that the effective action for p will also receive contributions from a Jacobian due to 
the change of variables from the eigenvalues to the density [|lT|. Since everything we have 
to say is independent of this term however we will just drop it here. 
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An equation for the eigenvalue distribution can be found by looking for stationary 
points of the effective action in Eq. (4). In performing infinitesimal variations of px one 
must be careful though not to change the normalization of p^ in Eq. (2). This can be 
achieved by taking variations of the form 

^P.(A) = (5) 

or equivalently 

5^{\) ~ NdX6p{\) ^ ^ 

where ip{X) vanishes at the end point of the support of p(A). The saddle point equation 
for (a translationally invariant) p follows by setting the first variation of the effective action 
S S 

— — — equal to zero. One obtains, denoting the logarithmic derivative of the Itzykson-Zuber 

6ip{X) 

determinant by F{X), 
where 

y'w = I d^'^- (8) 



A -A' 

This equation by itself does not determine p since F is unknown. The function -F(A) however 
satisfies a set of Schwinger- Dyson equations which were derived by Migdal in Ref . . Using 
these equations one finally ends up with an integral equation for Rel^'(A) whose imaginary 
part determines p: 

where 

RiX) = ^ReV'iX) + (10) 

In one dimension the gauge field can be gauged away and using U{^) = |m^$^ the action 
in Eq. (1) is just that of a free scalar field. In this case it is well known that p is semicircular 
in shape, i.e. 

P(A)=Po(A) = ^Jp-^ (11) 



and p = \/m^ — 4. Gross 0] realized that the semicircular form (0) for p actually solves 
Eqs. (9,10) in any number of dimensions. He showed that in D dimensions 



m^{p - 1) ± DJm^ - 4(2D - 1) 

= W^l (^2) 
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and also computed the free energy^: 

^ 1, I T , ,1 1 / T,, , , 

The following things are important to note: /i+ is a minimum of F and vanishes at m? = 2D 
in D < 1 implying that in this case the continuum limit can be taken. In D > 1 on the 
other hand, /i+ never vanishes and instead /i_ — >■ as — 2D~ . However the saddle point 
at /i = /i_ is a local maximum which, as we shall see, has important consequences for the 
spectrum of the theory. 

To summarize this section: For a quadratic potential the distribution of eigenvalues is 
semicircular in shape in any number of dimensions. In one or less dimensions a continuum 
limit can be constructed and we expect a physical spectrum since one is perturbing around a 
local minimum of the free energy. In more than one dimensions on the other hand one expects 
tachyons in the continuum limit. The spectrum of fluctuations around /i+ will describe the 
spectrum of mesons made out of adjoint matter in the strong coupling lattice theory. We 
will now go on to show all this explicitly. 



Migdal's Wave Equation 



In order to get at the spectrum of the theory in leading order in — one must work out the 
effective action describing the fluctuations around an extremum of the action (4). This has 



been done by Migdal [|Tl|]. Writing 

6p,{X)=p,{X) + 6p,{X) (14) 
with Sp as in Eq. (5) and PxiX) a solution of (7), one obtains 

S2[p] = -J2jdxJdX' ^-v{X, X')^|JM^|Jx+a,{X') + [j^^^^j^ + Da{X, X')]M^)M^') 

(15) 

In the case at hand here we have Px{X) = Po{X). Now using plane waves, ipxW = 4'{X)e^^^, 
one immediately obtains Migdal's wave equation for the particle spectrum: 



J dX'T^iX, AX A') = -J dX'lj^—j;^ + DaiX, A')]^(A'). (16) 
In Eqs. (15,16) we have introduced Migdal's notation for the derivatives of F{X), i.e. 



^It is interesting to note that the the free energy obtained by Kawamoto and Smit Q can be brought 

4 iTfi^ 

into the same form if one makes the curious substitution A = and M — , where A is the square of 

IJ,-' 8 

the meson field and AI is the bare quark mass. 
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and 



d SF{X) 



'1^ 



Also, in Eq. (16) we have defined Vt^ = I]^cos(P^). Migdal |Tl| has derived Integral 
equations for the functions rj and a. We have checked his derivation of these equations and 
agree up to an additional term in the equation for a. The origin of this term is explained in 
Appendix 1. Using the correct equations for t] and a from the Appendix we obtain, denoting 

— — ■?/;(A') = H{\') and integrating Eq. (16) over A, the following equation: 



/ d\ 



A'- An 



dX / dX'H{X')K{Xo, A) 



D)^^G{Xo,X') 1 



X' ' ^ G{Xo,X) A -A' 
-D I dX'H{X')G{Xo, X' 



(19) 



Several comments are in order. First, all the integrals in Eq. (19) are defined in the principal 
value sense. Second, Migdal's functions K{Xo, A) and G{Xo, A) are given by 



K{Xo,X) 



P(A) 



7rV(A) + (Ao-i?(A))^ 



and 



G(Ao,A) 



1 



Ao - R{X) 



Re exp 



+ 



dX' 



■ arctan ■ 



np{X') 



(20) 



(21) 



7r(A' - A) Ao - R{X') 

Third, the range of the integrals is over the support of p(A) which in our case is the interval 
[—^, +;^]- Fourth, the last term in Eq. (|19D is the one absent in Migdal's paper |1TT| . 



In principal the calculation of the spectrum of fluctuations around the Gross saddle point 
is now straight forward: Compute G(Ao,A) using it's definition Eq. (0), do the A-integral 
in ([19D and finally solve the resulting (singular) integral equation. However there are several 
subtleties in the calculation which we thought make it worthwhile to present the calculation 
of G{Xq, A) in some detail. This is done in Appendix 2. The interested reader will be able 
to check all the other integrals relevant here by using the method used in the Appendix. 
Suffice it to say here, that boundary conditions are important. For example, the branch of 
the 'arctan' in Eqs. (|T^,^T]) is chosen by a boundary condition at infinity 0]. Hence, in the 
calculation Aq must be taken outside the support of po- 



We have obtained the following expression for G{Xq, A): 

1 {RXo - A) - V^-p 



^(AcA) 



(^-f) (Ao-i?A)2 + p- 



(22) 



1 



where -R(A) = RX = ^^("^^ + f^i^ ~ -'-))A is linear in A for the semicircular solution. To 
derive our final result for the wave equation below, we have repeatedly used the important 
identity 



(23) 
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satisfied by n±. Finally, for completeness we also give the results of the A-integrals in Eq. 
(|19D. We obtained 

^^ K{\o,\) ^ fil{R-xy) 

and 



G(Ao,A') 1 u J{Rx- - I 

where we have introduced the new variables x = and y = ^A'. 



Calculation of Spectrum 

Using the results of the previous section we obtained the following integral equation deter- 
mining the spectrum of the theory to leading order in 

^ / . dy-^ = I dyh{y)- . J (26) 

where 



and 



-ly-x J-i ^x-RyY + i^{l-y^) 

A{x) = (l^£k£ - Dx (27) 

Bix) = ii^£te - DR. (28) 

To bring the wave equation into this form we have assumed that H{X') depends on A' only 
through y, i.e. H{X') = h{y). 

To determine the spectrum of the theory means to determine in Eq. (26). Due to the 
compactness of the integration interval we expect a discrete spectrum. To solve the equation 
we follow the time honored method of expansion in a complete, orthonormal set of functions. 
In particular, due to the range of the integration interval and the form of the integrand on 
the left hand side of the equation, Chebyshev polynomials are a natural set of functions to 
use here0. Hence we make the ansatz 

n V ^ ~ y 

Using this ansatz the integral on the left hand side of Eq. (26) can be done immediately 
(remembering though that |a;| > 1 in the calculation) and one obtains 



Q / dy = -7^^]2^c„ , „ ^ • (30) 

J-i y-x V va;2 - 1 



^We use the conventions of Rcf. [Oj for the definition of the Chebyshev polynomials. 
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On the right hand side we first write 

A{x) - B{x)y f{x) g{x) 



{x - Ryf + ^(1 - y2) y-y+ y - y_ 

where 



(31) 



y± = Rx± |Vx2 - 1. (32) 

In the last equation the identity Eq. (23) has been used. Now the integral on the right hand 
side of the wave equation (26) can also be done and one obtains after some algebra, 



'''\x-Ryr + ^{l-y 



-^EU[P-l) + ^](^+fr + [D-{D-l)]{R+^r}^^^^^^. (33) 

2 n ^ ^ Vx-1 

Now comparing Eqs. (30) and (33) we finally obtain for the eigenfunctions and eigenvalues 

hniy) = (34) 

and 

= ^-^Y^m + f -iR+ f )i + \r + iR+ f )-"]• (35) 

The last two equations constitute the main result of this paper. We will discuss it's sig- 
nificance in the next and final section. It is noted that the fluctuations themselves (at a 

1 T ( ^ 

given momentum P„) are of course given by 5p{y) = " which properly satisfies 

N ^yi — y^ 

J dySp{y) = for n = 1, 2, ■ ■ but 7^ for n = 0. Therefore ho{y) should be excluded 
from the eigenfunctions. 



Discussion and Conclusion 



As expected we have found a discrete spectrum with an infinite number of states. The 
quantum number n labels states according to their parity under y —y. The eigenfunctions 
hn{y) are even (odd) under this operation for n even (odd). Let us begin our considerations 
of the spectrum now by looking at the one dimensional case in some detail. To this end it 
is useful to define 

M = ln{R + ^) (36) 

and write Eq. (|35| ) as 

= D cosh(nM) - (D - 1) sinh(nM) (37) 
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In one dimension we obtain 



cos(P„) = cosh(n ■ arccosh(^-)) = T„(-^) (38) 

The interpretation of this result is straightforward. As we mentioned before, in one dimension 
the gauge fields play no role and the theory is free with propagator 

1 1 
N — 2 cos(P) 

Clearly, the 'pole' of the propagator is at cos(P) = ^ which agrees with (38) in the case 
n = 1. In position space this corresponds to an exponential falloff of the propagator with 
mass M. The states with n > 1 correspond to 'mesons' made out of n $'s. The x-space 
propagator for these objects decays exponentially with mass n-M which in momentum space 
gives Eq. (38). 

In one dimension /i = ^/mF^^ vanishes as m? — > 2+ and we can take the continuum 
limit. Writing = 2 + mpa^ and P = Pga one obtains 

- ^El = v?. ml (40) 

Hence in the continuum limit we get a linearly rising spectrum of states of mass nmQ. Note, 
that for mo = we get a massless boson in one dimension as opposed to the massless boson 
in 1 + 1 (one spacelike, one timelike) obtained in the double scaling limit of the continuum 
model ll^. In less than one dimension still vanishes at th? = 2D and one obtains expanding 
Eq. (^) at = 2D + m^a^ 

El = n-ml (41) 

in the continuum limit. Contrary to the case in one dimension, the particle masses are 
proportional to ^/n. 

Let us now go on and discuss the case D > 1. As we pointed out before the stable solution 
is n = in Eq. (12) which is always greater than zero. In this case the interpretation of 
the spectrum in Eq. ( pTf ) in four dimensions is that of the mesons in strong coupling lattice 
QCD with adjoint matter. Note that asymptotically in n, = (-R + /^/2)". The same 
formula holds in the limit D ^ oo for all n. It is interesting to see what happens if instead 
we use fi = fi-. Setting = 2D — m^a? in order to approach the critical point from below, 
one gets 

-P^^^El = ~n-ml (42) 

in the continuum limit. Although infinitely many states appear in the continuum limit, they 
are all tachyonic as expected. The spectrum in D > 1 could be improved if an upside-down 
quadratic potential[0 is used for t/($) in Eq. (|I]). The continuum limit in this case is 
equivalent to setting = —2D — rn^a? in Eq. (0), which gives 

Ey = 2D ((-1)" - 1) + (-1)" ■ n ■ my. (43) 

For n even, we obtain En = n ■ m^ > though En = —oo for n odd. 
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What do the above results tell us? In one dimension we got the expected result. In more 
than one dimension we find an infinite tower of states which certainly suggests a 'stringy' 
interpretation of the spectrum. Makeenko has written down Schwinger-Dyson equations 
for the functions 

G^C^y) =< ltr($(x)f/(a,)$(y)f/(C^,.)) > (44) 

where U{Cxy) denotes a string of links along the path C^y For the quadratic potential 
Makeenko has obtained the solution for G{Cxy) and it would be interesting to see whether 
by summing his result over all paths our spectrum can be reproduced If it is an 

interpretation of our result in terms of the excitations of a string of fiux seems natural. 
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Appendix 1. Derivation of the extra term in eq.([T9|) 



In this Appendix we show how the last term in eq. (^) is derived. In the Appendix of Ref . 



ni Migdal has derived the following equation (78) 



p(/x)^A(0^/i) /'^-A ^^'^) (^('^'''^) + j^^'?</>'-S'^'^^) 

z-^i -^^^^^i ^^^^-0 (A-/?(0))2 + vrV(0) ^ ' 

where 

Setting z — > cxD in the above formula, we find 



From eq.(59) in Ref. |jrT 



/ df,pMGxifi) = J du^f^ , (4^ 



we obtain 

df,p{fi)g^{(j)', /i) = ~^Gx{<l)') , (49) 



^In one dimension one just obtains the square of the x-space propagator for a free particle by this 
procedure, i.e one only obtains the lightest 'meson'. 
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which differs from eq.(77) in Ref. [O]. Inserting 
a(0',0) 







into (^7D, we finally find equation for 

(50) 



(A - + 7rV(0) 



Appendix 2. Calculation of ^(Ao, A) 

The definition of G{\o, A) is 

1 [ r 'rrn(\'\ 1 

(51) 



rfA' 7rp(A') 
exp[/ — — — - arctan ■ 



7r(A'- A) Ao-i?(A')^ 



where -R(A) = -RA for the semicircular solution 7rp(A) = \J jJ' — We first calculate the 

following integral 

J = I +^^^ arctan " ^"^"'^ (52) 
Integration by parts together with the obvious rescaling of variables gives 

J=^l£/'^l„(,-.-)|^ + ^] (53) 

TT 7-1 VI - r y-y+ y- y- 



where x = ^Aq, z = ^A, y = ^X', f± = ±— — , and y± are given in Eq. 

y+-y~ 

the text. Using the following properties for Chebyshev polynomials T„(?/) 



m 



oo /j. \n 

ln{y -z) = -To{y) ln(2t,) - 2 ^ T^{y) ^ (54) 



1 _ 2t± 

72 



y-y± 1 - 



To(i/) + 2^r„(t±)") (55) 



± n=l 



'^;i^r.(y).r,„(j,)=^„„x{-/2 (56) 



where = 2; — iVl — z'^ and t± = (-R =F f )(3; ~ Vx'^ — 1), we perform the integral in Eq. 
)3D and obtain 

1 - 



J=-\n[ ^ 1 . (57) 



From the above result G becomes 



1 „ . 1 - t,t+ 



which finally gives eq.(P^ in the text. 
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